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ABSTRACT 

We present a study of H II regions in M51 using the Hubble Space Telescope 
ACS images taken as part of the Hubble Heritage Program. We have catalogued 
about 19,600 H II regions in M51 with Ha luminosity in the range of L = 
10 35 ' 5 erg s~ x to 10 39 erg s _1 . The Ha luminosity function of H II regions 
(H II LF) in M51 is well represented by a double power law with its index a = 
— 2.25±0.02 for the bright part and a = — 1.42±0.01 for the faint part, separated 
at a break point L = 10 37 1 erg s -1 . This break was not found in previous studies 
of M51 H II regions. Comparison with simulated H II LFs suggests that this 
break is caused by the transition of H II region ionizing sources, from low mass 
clusters (with ~ 10 3 M , including several OB stars) to more massive clusters 
(including several tens of OB stars). The H II LFs with L < 10 37 1 erg s -1 are 
found to have different slopes for different parts in M51: the H II LF for the 
interarm region is steeper than those for the arm and the nuclear regions. This 
observed difference in H II LFs can be explained by evolutionary effects that H II 
regions in the interarm region are relatively older than those in the other parts 
of M51. 

Subject headings: galaxies: individual (M51; NGC 5194; NGC 5195) — galaxies: 
spiral — galaxies: ISM — H II regions 



1. Introduction 



H II regions are an excellent tracer of recent star formation since the radiation from H II 
regions carries a signature of young OB stars. The total flux of the hydrogen recombination 
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line (e.g., Ha) coming from an H II region is proportional to the total Lyman continuum 
emission rate of the ionizing stars, if we assume that all ionizing photons are locally absorbed. 
Thus, the Ha luminosity of the H II region (L) is an indicator of the Lyman continuum 
emission from ionizing sources. 

In general, H II re gions are classified into three groups, depending on their observed 
Ha lu minosity and size flHodgelll969l . Il974l iKennicutt etaDll989l : iHodge et al.lll989l : iFranco et al 



20041 ). (a) Classical H II regions with Ha luminosity fainter than L « 10 37 erg s 1 : With 



their sizes up to several parsecs, these H II regions are typical H II regions ionized by several 
OB stars. A representative example of this class is M42, the Orion nebula, (b) Giant H II 
regions with Ha luminosity of about L = 10 37 erg s _1 - 10 39 erg s -1 : These H II regions are 
ionized by a few OB associations or massive star clusters and usually smaller than 100 pc. 
Typical examples in our Galaxy of this class are W49, NGC 3603, and the Carina nebula, 
(c) Super giant H II regions with Ha luminosity brighter th an L ^ 10 39 erg s~* : These 



may be ionized by multiple star clusters or super star clusters (IWeidner et al.l 120101 ) . These 
H II regions have no known analogue in the Galaxy, and they are mostly found in late-type 
galaxies or interacting systems. Two examples of this kind are 30 Doradus in the Large 
Magellanic Cloud (LMC) and NGC 604 in M33. 

The Ha luminosity function of H II regions (H II LF) in a galaxy provides a very useful 
information on the formation and early evolution of stars and star clusters, and has been 
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Kennicutt et al. 1989; Rozas et al. 



Bradley et al.l 120061 ). It is usually represented 



by a power law of the following form: 



N(L)dL = AL a dL 

where N(L)dL is the number of H II re gions with Ha luminosi ty between L and L + dL, a is a 
power law index, and A is a constant. IKennicutt et al.l (119891 ) found that the Ha luminosity 
functions of bright H II regions in 30 nearby spiral and irregular galaxies are represented 
approximately by a power law with a = —2.0 ± 0.5. For some galaxies, H II LFs are known 
to have a break at L fa 10 38,9 erg s~\ being steeper in the bright part than in the faint part 



( IKennicutt et al.l Il989t iRandl Il992l ; iRozas et al. 



19961 : lKnapenlll998f ). Th e H II LFs with 



this break is called 'Type II', and this break is called Stromgren luminosity (IKennicutt et al. 



19891 ; iBradley et al.ll2006l ). The H II LF slope also depends on the reg ions in a spiral galaxy : 



the steeper H II LF for the interarm region than fo r the arm region ( IKennicutt et al.l 11989 
Randlll992l : iThilker et alfaoot IScoville et aDboOlh . 



Several scenarios have been suggested to explain the cause of the H II L F slope change at 
10 38 9 erg s -1 : different molecular gas cloud mass spectra ( IRandlll992[ ). the evolution of 
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H II regions and the variation in the number o f ionizing; stars ( Ivon Hippel fe Bothun 111990 



McKee fe Williamsl Il997t lOev fe C larke 



density-bounded H II regions (IBeckman et al 



1998 



caused by low resolution observations fjPleuss et al 



the transition from ionization-bounded to 



2000 ), an d the blending of small H II regions 



20001 ). All these scenarios, except for one 



involving the observational effect, basically suggest that the H II LF shape should depend 
on the physical conditions of the corresponding regions in a galaxy. 



Most studies on the H 



fract i on of faint H II region s fcennicutt fe Hodge!ll98(| 



200CH : iThilker et all 120021 : Bradley et al 



I LFs are based on ground-based images, missing a significant 



2006 



Kennicutt et al.l 1989; Bcckman et al. 



Helmboldt etaD 12009^ except for a sm all 



number of the galax ies in the Local Group fcennicutt fe Hodgelll986l : iHodge et alill989l lbl: 



Hodge fe Lee 1 1 1 9 9 Oh . Therefore, the nature of the faint H II regions in galaxies is not well 
known. There are a few studies based on Hubble Sp ace Telescope (HST) images, but they 



covered only a small fraction of the target galaxies (IPleuss et al.ll2000l ; IScoville et al.l 12001 



Buckalew fe Kobulnickyj|2006l ). We started a project to study the H II regions in M51, using 
the high resolution images taken with the HST Advance Camera for Survey (ACS). 

M51 is a system of interacting galaxies, including a grand-design spiral galaxy NGC 
5194 (Sbc) and a b arred lenticular galax y NGC 5195 (SB0), located in a relatively close 



dista nce (9.9 Mpc, iTikhonov et al.l 120091 ) . Its almost face-on orientation (i = 20°, iTully 
1974 ) enables us to observe its structure in detail with minimal obscuration by interstel- 
lar dust. Spiral structures in M51 have been used as a s trong constr a int to model the 



origin of spiral arms and the effect of galaxy intera ction (IDobbs et al. d2010h and refer- 



ence t herein). M51 is ab undant in interstellar media (ISchuster et al. 



2007 



Hitschfeld et al. 



20091 ; lEgusa et al.l 120101 ) and i s well known for ac tive star formation that has sustained 



at least longer than 100 Myr (ICalzetti et al.l l2005f ) . This is consistent with the existence 



of m a ny young and intermediate age star clusters (IHwang fe Led 120081 ; IScheepmaker et al. 



2009; 


Hwang & Lee 


2010: 


Kaleida & Scowen 


201 



1989; 


Rand 


1992; 


Petit et al. 


1996: 


Thilker et al. 


200(1 



201o[) as well as numerous H II regions along 



the spiral arms of M51 (|Carranza et all Il969l: Ivan der Hulst et aP Il988l; iKennicutt et al.l 



However, most of the existing studies on M51 H II regions are based on the ground-based 
observations. In a study of 616 H II regions with L > 10 36,4 erg s _1 in M51, iRandl (119921 ) 
reported that the H II LF has an inverse slope in the L < 10 37 6 erg s _1 , and that the H II 
LF slope changes at L « 10 39 erg s _1 . He also found that the H II LF slope for the interarm 
region (a = —2 .05 ± 0.15) is steeper than that for the arm region (a = —1.48 ± 0.07). 



Petit et al.l (119961 ) estimated the Ha luminosities and sizes of 478 H II regions in M51 using 
the Ha photograph taken with the Special Astronomical Observatory 6 m telescope, and the 



Ha luminosity ranges from L = 10 36 ' 4 erg s 1 to 10 39 - e erg s 1 . They reported a nearly flat 
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slope below 10 37,6 erg s 1 , in contrast to the result of iRandl (119921 ). and the H II region size 
distribution is well fitted by an exponential function. 



Later, iThilker et al.l ( 120001 ) found about 1,200 H II regions from Ha image data covering 
a central 6.7' x 6.7' field of NGC 5194 using HIIPHOT that was developed for the detection 
and photometry of H II regions. They showed that there is a change in the slope of the 
H II LF at L — IP 38,8 erg s _1 , consistent with the result of IRandl (1l992l ). On the other 



hand, IScoville et al.l (120011 ) detected about 1,400 H II regions using the HST Wide Field 
Camera 2 (WFPC2) observation covering the central 4.5' x 4' field of NGC 5194. The H II 
LF turns out to be steeper for the interarm region (a = —1.95 ± 0.05) than for the arm 
region (a = —1.72 ± 0.03). Howe yer, they cou ld no t find any change in the H II LF slope at 
L = 10 38 - 8 erg s _1 , in contrast to IRandl Jl992h and IThilker et all (boooh . 



These previous stu dies o n the H II regions in M51 h ave so me limitation s. The ground- 
based images used by IRandl (jl992j), iPetit et al.l (119961 ). and IThilker et al.l (120001 ) cannot 
resolve blended H II regions into small or faint H II regions due to its low spatial resolution 
of about 1.8" or 100 pc at the distance of M51. Also, it is difficult t o distinguish b l ended 
low-luminosity H II regions from diffuse ionized gas. On the other hand, IScoville et al.l (120011 ) 
used HST image data with 0.1" — 0.2" resolution 4.8 — 9.6 pc at the distance of M51), 
which enabled to resolve blended H II regions and to distinguish faint H II regions from 
diffuse ionized gas. However, their data cover only the central part of NGC 5194. 

Recently, M51 was observed with the HST ACS in Ha band as well as other continuum 
filters as par t of the Hubble Herit age program, and the data were released to the astronomy 
community (IMutchler et al.ll2005l ). These deep and wide field Ha imaging data provide an 
excellent opportunity to study the global properties of H II regions in M51. In this study, 
we present a result of H II region survey over the field covering most of NGC 5194 and NGC 
5195 regions using these data. We make a catalog of resolved H II regions, derive the H II 
LF, and investigate the properties of H II LF in different parts of M51. 



During this study, iGutierrez &: Beckmanl (120101 ) presented a study of over 2000 H II 
regions detected in the same HST ACS image data and reported the dependence of the 
mean luminosity weig hted electron density o f the H II regions on both H II region size and 
galactocentric radius. IGutierrez et al.l (120111 ) presented a catalog of the 2,659 H II regions 
in M51, and showed their Ha luminosity function and size distribution. 

This paper is composed as follows. §2 describes the data, and §3 introduces the H II 
region detection and photometry procedures. §4 presents a catalog of M51 H II regions. We 
analyze the properties of the H II regions, including H II LF and size distribution in §5, and 
discuss the primary results in §6. Finally, a summary and conclusion is given in §7. 
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Data 



M51 was observed with th e HST ACS in January 2005 through the Hubble Heritage 
program ( IMutchler et al.l 120051 ). The data set is composed of images in four filters ^43514^ 
(B), F555W (V), F8UW (I), and F658iV (Ha) with the accumulated exposure times of 
2720, 1360, 1360, and 2720 seconds, respectively. All the necessary data reduction processes 
were conducted by the Hubble Heritage Team, including the multi-drizzling for image combi- 
nation before the public release of the data. One pixel of HST ACS mosaic data corresponds 
to 0.05" after the multi-drizzling process. The FWHM of a point source is about 0.1". The 
size of the field of view is 7.2' x 10.2'. The spatial coverage is large enough to include the 
entire disk of NGC 5194 as well as its companion NGC 5195. Throughout this study, we 
adopt a distance to M 51 of 9.9 Mpc derived using the tip of the red giant branch method 
( ITikhonov et al.l 120091 ) . At this distance, one arcsecond corresponds to a linear scale of 48 
pc. 



3. H II Region Detection and Photometry 

Detection and photometry of H II regions involve following steps : (1) continuum sub- 
traction from the Ha images, (2) source detection and flux measurement on the continuum- 
subtracted images, and (3) correction for the [N II] line contamination and the extinction to 
derive net Ha fluxes of the H II regions. 

The effective bandwidth of the F658N filter is 74.8 A, so that the images obtained with 
this filter include not only the Ha line emission but also the [N II] 6548A, 6583A line emission 
and the continuum emission. We need to subtract the continuum from the F658N images to 
find the Ha line emitting objects, and to measure their Ha emission line fluxes. We used the 
average of F555W and F814W images to make a continuum image. We measured the flux of 
22 non-saturated foreground stars in the F658N image and the combined continuum image 
and derived a scaling factor 0.0878 from the average ratio of these intensities. Then, the 
combined continuum image was multiplied by the derived scaling factor, and was subtracted 
from the F658iV filter image to make a continuum-free Ha image. Figure [T] displays the 
continuum-subtracted Ha image of M51, showing bright, discrete H II regions outlining the 
spiral arms of NGC 5194, and the significant amount of diffuse Ha emission in the region of 
NGC 5195. 



We detected the H II regions and determined their fluxes and sizes in the continuum- 
subtracted im age using HIIPHOT that was developed for the detection and photometry of 
H II regions by lThilker et al.l (120001 ). HIIPHOT determines the boundaries of individual H II 
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regions using the parameter that dictates the terminating gradient of the surface brightness 
profile. We adopted a value for the terminal surface brightness gradient, 10 EM pc _1 where 
EM is Ha emission measure in the unit of pc cm -6 . One EM corresponds to a surface 
brightness of 2 x 10~ 18 erg cm -2 s _1 arcsec -2 . The local background level for each H II 
region is obtained using a surface fit to pixels located within an annulus with 4" width 
outside the boundary. We used S/N > 12 as a detection limit. With this condition we 
found 19,598 H II regions. We calibrated the instrumental Ha fluxes (fu a ) of the H II 
regions using the photometric zero point (PHOTFLAM) given in the header of F658N image, 
1.999 x 1CT 18 erg cm -2 A -1 . Then we derived their Ha luminosity using L = And 2 x in a 
where d is a distance to M51, 9.9 Mpc. 

We applied two corrections before obtaining the final photometry: [N II] contamination 
and the Galactic foreground extinction. First, the Ha flux in the continuum-subtracted 
image includes contributions from the two satellite [N II] lines. At zero redshift, the [N II] 
lines are located on either side of the 6563A Ha line, one at 6548A and the other at 6583A. To 
measure the net Ha flux, we consider the redshift (~ 0.002475) of M51, transmission values of 
the HST ACS F658iV filter, and me an [N II] /Ha ra t io (~ 0.5) obtained from spectroscopic 
data for ten H II regions in M51 by iBresolin et all (120041 ) . We used Ha / (Ha + [N II]) 
~ 0.67 to derive the net Ha flu x. Second, we applie d the foreground reddening correction to 
the Ha flux using the values in lSchlegel et al.l (119981 ) : E(B-V)=0.035. For the corresponding 
extinction (Ay = 0.115 mag, Ah q = 0.092 mag), the measured Ha flux is increased by a 
factor of 1.09. Internal extinctions for individual H II regions are not known. We adopted 
the mean extinction value Ay ~ 3.1 mag derived from Ha/Paa flux ratio for 209 H II regions 
in the central region of M51 by lScoville et al.l (120011 ) . This value corresponds to a change in 
Ha luminosity by a factor of ~ 10. 

There are some sets o f Ha photometry d ata for the M 5 i H II regions avail a ble in 
the literature: 



Rand 



and Gutierrez et al 



19921 ). IPetit et al.l (Il996f l. iThi 



ker et al. 



Randl rtl992h 



jgOOol). IScoville et al.l ffipOlh 



Petit et al 



(Il996h 



and 



20111 ) . Since the studies by 
Thilker et al.l (l2000f ) are based on the ground-based images, it is difficult to compare di- 
rectly our photometric measurement s with theirs. There fore, we compared our res ult of 
H II region photometry with those of IScoville et al.l (120011 ) and iGutierrez et al.l ( 1201 ll ) . We 
selected 54 isolated , compact, and bright H II regions common between this study and 
Scoville et al.l (120011 ) by visual inspection, displaying the comparison of the measured fluxes 
in Figure ET a) . Figu re BTa) shows that our photometry is in good a greement with that of 
Scoville et al.l ( 1200 ll ). The mean difference (log L [this study] - log L [Scoville et al.l ( 120011 )]) 
is 0.03 ± 0.12. For Figure |^b), we selected 164 i s olated , compact, and bright H II regions 
common between this study and IGutierrez et al.l (120111 ) by visual i nspection. Figu re Mb) 
shows that our photometry is also in excellent agreement with that of IGutierrez et al.l ( 1201 ll ). 
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The mean difference (log L [this study] - log L Gutierrez et al.l (]201lh ]) is 0.01 ± 0.07 



4. H II Region Catalog 

Measured properties of H II regions in a galaxy are dependent on both the resolution of 
images and the procedures used to define the boundaries of H II regions. For example, higher 
resolution images can resolve some large H II regions into multiple components, reducing the 
number of the detected large H II regions. In addition, the criteria adopted to separate the 
H II regions from the background are critical. In this study, we present two H II region 
catalogs: for the original and group samples. The original sample includes H II regions 
detected and provided from HIIPHOT, and the group sample includes H II regions produced 
by combining smaller ones in close proximity into one large one. 

The original sample catalog includes about 19,600 H II regions, containing resolved 
sources in blended H II regions as well as isolated H II regions. This catalog is useful for the 
study of the relation between classical H II regions and star clusters since most H II regions 
in the original sample are classified as classical H II regions, and only ten percent of them are 
giant H II regions. There is no super giant H II regions in the original sample. In Table CD, we 
list positions, Ha luminosities, and effective diameters (derived from D = 2 x (aiea/n) 1 ^ 2 ). 
Table [T] lists a sample of H II regions for reader's guide and the full catalog will be available 
electronically from the Online Journal. 

To prepare the group sample, we found as sociations of H II regions in the original sample 



using the Friends-of- Friends algorithm (FOF. lDavis et al.lll985l ). Every source in groups has 



a friend source within a distance smaller than some specified linking length. If we increase 
the linking length, the number of members in individual groups increases, decreasing the 
number of groups. As a test, we made 20 group sample sets using the linking length ranging 
from 0.1 to 2.0 mean separation length with a step of 0.1. Through this test, we found 
that the linking length of 1.0 mean separation is the most suitable in reproducing large H II 
regions found in ground-based images. 

We call these associations of H II regions as grouped H II regions. The total number 
of the grouped H II regions is 2,294. We inspected these grouped H II regions on the 
continuum-subtracted Ha image, revised the membership of 209 H II regions. The final 
number of grouped H II regions is 2,296, and the number of the isolated H II regions which 
do not belong to any group is 4,919. Thus, the sum of these two kinds of H II regions is 
7,215. The position of a grouped H II region is derived from average values of the positions 
of all individual H II regions in a group. The area and luminosity of a grouped H II region 
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are derived from summing the values of all individual H II regions in a group. Most of the 
2,296 grouped H II regions are giant H II regions with L > 10 37 erg s _1 , and 18 of them are 
super giant H II region with L > 10 39 erg s _1 . Table [2] lists positions, Ha luminosities, and 
effective diameters of H II regions in the group sample. The last column in Table [T] marks 
the group ID numbers. 



We compa red positions of 7 ,216 H II regions in the group sample with those of 2,659 
H II regions in I Gutierrez et al.l (120 111 ). 2,402 HII regions in the Gutierrez et al's catalog 
were found in our catalog, showing that about 90 per cent of the Gut i errez et al's catalog 
is included in our catalog. Among the H II regions in iGutierrez et al.l ( 1201 ll ) not matched 
with our H II regions, 63 H II regions were included as part of the grouped H II regions in 
the group sample, and 95 H II regions were splitted into multiple H II regions. The genuine 
membership of these H II regions is very difficult to determine. 20 H II regions are located 
outside the field used in our study. There are 79 H II regions not found in our survey. We 
inspected carefully the continuum-subtracted Ha image of these H II regions, and found (a) 
that there is no recognizable Ha emission for 53 H II regions, and (b) that 26 H II regions 
are located in the faint diffuse emission feature, but with lower S/N than used in our survey. 



5. Results 



5.1. Ha Luminosity Function of H II Regions 



The total number of M51 H II regions listed in the original sample is about 19,600. 
The Ha luminosity of these H II regions ranges from L = 10 35 5 erg s" 1 to 10 39 - erg s _1 . 
This detectio n limit is similar to the deep survey of H II re gi ons in the dwarf gal axies of the 
Local Group fcennicutt & Hodgelll9~8~ol : iHodee et aHll989l Jbl: iHodge fc Lee Ill99ph . The total 
Ha luminosity of all these H II regions is derived to be L = 10 4L2 erg s -1 . There are 3,653 
H II regions with L > 10 37 erg s _1 , and their total luminosity corresponds to 70 percent of 
the total Ha luminosity of M51. There are only 160 H II regions with L > 10 38 erg s _1 . 

We derived the H II LF for the original sample, as shown in Figure |3j The lower limit 
of Ha luminosities is determined by the observational threshold on the surface brightness. 
The H II LF shows a maximum at L fa 10 36 1 erg s _1 and declines rapidly as the luminosity 
decreases. This value is similar to the o ne for M33 H II regi ons, L fa 10 35 - 9 erg s _1 , found from 
the completeness-corrected H II LF bvlWvder et al.l (|1997l). and those for N GC 6822 and the 



Magellanic Clouds, L fa 10 36 erg s 1 ( IHodge et al.lll989t IWilcots et al.lll99ll ). Correcting this 



value for the adopted mean extinction value for M51 H II regions leads to ~ 10 37 1 erg s 
This value is similar to that for M42, fa 10 36 9 erg s~ 1 (jScoville et al.ll200ll ). Considering that 
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the H II LF for M33 shows a flattening for L = 10 35 erg s' 1 - 10 36 erg s" 1 ( Hodge et al.lll999f ). 
the rapid decline for L < 10 36 erg s _1 in our H II LF appears to be due to incompleteness in 
our catalog. 

The H II LF appears to be fitted by a double power law with a break at L ~ 10 37 ' 1 erg s -1 . 
We used a double power law function to fit the H II LF: 



N(L)dL = AL ai dL a\ : Bright power law index for L > Lb 



and 



N{L)dL = A'L a2 dL 



where is the break point luminosity, and A' = AL^ 2 



«2 : Faint power law index for L < 

For the bright H II regions 



N erg s 1 ), we obtained a power law index, a = —2.25 ± 0.02. 



(10 37 - 1 erg s- 1 < L < 10 
For the faint part (10 36 ' 1 erg s _1 < L < 10 37 ' 1 erg s -1 ), the H II LF becomes flatter than for 
the bright part, having a = —1.42 ± 0.01. 

A distinguishable point with this H II LF is a break at L = 10 37 1 erg s -1 . This break has 
not been observed in previous studies of the H II r egions in M 5 1, because the data used in pre 
vious studies were not de ep enough to detect it ( ]Randlll992l ; iPetit et al.lll996l ; iThilker et al. 
2000l ; IScoville et al.ll200ll ). Howeve r, the H II LF break at sim ilar lumino sity is known to ex - 
ist for some nearby galaxies: M31 (IWalterbos fc Braun Ill992l ) and M33 (IHodge et al.lll999l ). 
This break may be caused by the transition of H II region ionizing sources, from low mass 
cluster s (including several OB stars) to more massive clusters (including sev eral tens of OB 
stars) fcennicutt et all Il989l iMcKee fc Williamsl Il997l : bev fc Clarkel Il998h . This will be 
discussed in detail in the next session. 



We compared the H II LF for the original sample (solid line) with that in IThilker et al. 



( 120001 ) (dashed line), as shown in Figure HJa). The Ha luminosity of H II regions for the origi- 
nal sample range s from L = 1Q 3 5 5 erg s _1 to 10 39 erg s , while it ranges L = 10 36 ' 4 erg s _1 to 
Thilker et al.l ( 120001 ) . The minimum lu minosity of the H II regions in the 



10 



39.4 



erg s 1 m 



original sample is about ten times fainter than that in Thilk er et al. feoool) . While this 
study shows a break at L — 10 37 1 erg s _1 , the study by IThilker et al. (hoool ) shows a break 



at L = lCr 8,8 erg s , much brighter than that in this study. The power law indices in this 
study ar e a = —2.28 ±0-02 f or the bright part and a = —1.38 ± 0.01 for the faint part, and 
those in IThilker et~aD fcoooh are a = -3.41 ± 0.08 and a = -1.61 ±0.03. 



Figure ID^b) shows the H II LF for the original sample (solid line) and that in lScoville et al 



( 120011 ) (dashed line). The spatial coverage of Ha images used in this study is large enough to 
include most o f the entire disk of NGC 5194 as well as its companion NGC 5195. However, 
Scoville et al.l (120011 ) obtained the H II LF from Ha images covering only the central 4.5' x 4' 
field of NGC 5194. For comparison, we plotted two H II LFs derived for the same area. The 
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minimum and maximum luminosities of the H II regions in two samples are similar, but there 



is a difference in the numbers of H II regions with L < 10 37 ' 1 erg s 1 . The ACS i mages used in 



this st udy have longer exposure time than that of the WFP C2 images used in IScoville et al. 



(120011 ) , so that we found many more faint H II regions than IScoville et al.l ( 120011 ) . Both H II 
LFs are fitted by a power law for the bright part (10 37,0 erg s _1 < L < 10 38,6 erg s _1 ). The 
H II LF for the original sample is s lightly steeper with a = —2.24 ± 0.03 than the H II LF 
with a = -1.91 ± 0.04 reported bv lScoville et al.l (boOlh . 



Figure HJc) shows the H II LF for the original sample (solid line) and that in lGutierrez et al. 



(120 111 ) (dashed line). The Ha luminosity of the H II regions for the original sample ranges 
from L = lO 35,5 erg s" 1 to 10 39 - erg s _1 , while it ranges L = 10 35 8 erg s -1 to 10 40 erg s _1 in 
(2011). While this study shows a break at L = 10 37 1 erg s" 1 , the study by 



Gutierrez et al 



Gutierrez et al 



(120 111 ) shows a break at L 
study. The power law indices in this study are a 
a = - 
and a 



10 38 - 7 erg s~\ 



much brighter than that in this 
2.21 ± . 02 fo r the bright part and 
1.41 ±0.02 for the faint part, and those in lGutierrez et al.l (120 111 ) are a = -2.29±0.17 
= — 1.67±0.03. This large discrepancy between the two studies is pri marily due to two 



factor s: (a) we considered the small local peaks for the large H II regions in iGutierrez et al. 



( 1201 it ) as separated H II r e gions in this study, (b) we found many isolated faint H II regions, 
missed in IGutierrez et al.l (120111 ) . 

In Figure |5^a) we compared the H II LF for the group sample (solid line) with that for 
the original sample (dashed line). Although both H II LFs are fitted by a double power law, 
the location of a break and the power law indices are different. The locations of the break 
are L = 10 37 1 erg s _1 and 10 38,8 erg s _1 , respectively. The power law indices for the bright 
part are a = —2.21 ±0.02 and — 2.60±0.21, and those for the faint part are a = — 1.41±0.02 
and —1.69 ± 0.01 for the original and group samples, respectively. 

Figure |^b) shows the H II LFs for t he grouped H II reg ions (dotted line), the group 
sample (solid line) and the H II regions in lThilker et al.l (l2000f ) (dashed line). The locations 
of the break and the power law indices for the latter two samples a re similar. The locations 
of the break for both the group sample and iThilker et al.l (l2000f ) are L = 10 38 8 erg s _1 . 
The power law indices for the group sample are a = —2.60 ± 0.2 1 for the bright part and 
a = -1.69±0.01 for the faint part, and those in lThilker et aD pXHJ) are a = -3.41±0.08 and 
a = —1.61 ±0.03. In conclusion, the H II L F break seen at L = 1Q 38 8 erg s" 1 in the previous 
studies based on the ground-based images (jRandlll992l ; iPetit et al.lll996l ; IThilker et al.ll2000l ) 
can be explained by the blending of multiple H II regions. 



In Figure El^c) we compared the H II LFs for the group sample (s o 



H II regions (dotted line), and the H II regions in IGutierrez et al.l (j201l[ ) (dashed line) 



id lin e), the grouped 



The H II LFs for the latter two samples show nearly the same shape. The locations of the 
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break and the power law indices for the group sample and in IGutierrez et al.l (1201 



arc 



similar. The location of the break for the group sample and that in IGutierrez et al.l (1201 if ) 
are L = 10 38 - 8 erg s _1 and 10 38 ' 7 erg s _1 , respectively. The power law indices for the group 
sample are a = —2.60 ± 0.21 for th e bright part and a = —1.69 ± 0.01 for the faint part, 
and those in IGutierrez etaD f)201lh are a = -2.29 ± 0.17 and a = -1.67 ± 0.03. Both 
studies show a good agreement at L > 10 37 ' 4 erg s _1 , however, there is a large difference 
in the number distribution of H II regions at L < 10 37 4 erg s -1 . It is primarily due to the 
existence of many iso lated faint H II regions that were found in our survey, but missed in 



Gutierrez et al.l ( 1201 if ). 



5.2. Spatial Variation of Ha Luminosity Function for H II Regions 

We investigated the spatial variation of the H II LF for the original sample. We separated 
H II regions in the original sample into four groups, according to their positions in the galaxy: 
arm, interarm, and nuclear regions of NGC 5194, and NGC 5195. The positions of the H II 
regions are listed in Table HJ There are 12,245 H II regions in the arm region, 4,422 H II 
regions in the nuclear region, and 2,639 H II regions in the interarm region of NGC 5194. 
We found only 292 H II regions in the region of NGC 5195. The total luminosity of the 
arm H II regions is derived to be L = 10 411 erg s _1 and that of the nuclear H II regions 
is L = 10 40 ' 6 erg s _1 , while that of the interarm H II regions is only L = 10 40 1 erg s _1 , 8 
percent of the total Ha luminosity of all NGC 5194 H II regions. 

Figure |6]Ja) shows the H II LFs for the arm (solid line), interarm (dashed line), and 
nuclear (dotted line) regions of NGC 5194. The minimum luminosities (L m 10 35 5 ) in the 
H II LFs are nearly the same in three regions, but the maximum luminosities are different 
depending on the position. The maximum luminosities for the arm and nuclear regions are 
L = 10 39 - erg s -1 and L = 10 38 8 erg s _1 , respectively. However, there is no H II region 
brighter than L = 10 38 ' 2 erg s _1 for the interarm region. 

These H II LFs are fitted by a double power law with a break at L = 10 37 erg s _1 . The 
power law indices for the bright part are a = —1.92 ± 0.03, —2.01 ± 0.09, and —2.08 ± 0.05 
for the arm, interarm, and nuclear regions, respectively. On the other hand, the power law 
indices for the faint part are a = —1.33 ± 0.01, —1.53 ± 0.05, and —1.33 ± 0.03. Although 
the H II LFs for the bright part (L > 10 37 erg s _1 ) have similar power law indices in three 
regions, the H II LF for the faint part (L < 10 erg s -1 ) is steeper in the interarm region 
than in the arm and nuclear regions. 

Figure El^b) and (c) shows the spatial variation of the H II LFs for the grouped H II 
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regions and the group sample, respectively. The H II LFs for the grouped H II regions in 
Figure Mjo) are truncated at L — 10 37 erg s _1 . We used a single power law to fit the H II 
LFs for the range from L = 10 37,0 erg s _1 to L — 10 38 ' 4 erg s _1 . The power law indices 
are a = —1.57 ± 0.04, —1.74 ± 0.08, and —1.45 ± 0.06 for the arm, interarm, and nuclear 
regions, respectively. The H II LF is steeper in the interarm region than in the arm and 
nuclear regions. On the other hand, the H II LFs for the group sample in Figure [6]^c) are 
truncated near L = 10 36 erg s -1 . The power law indices are a = —1.78 ± 0.04, —1.95 ± 0.08, 
and —1.58 ± 0.06 for the arm, interarm, and nuclear regions, respectively. The H II LF is 
also steeper in the interarm region than in the arm and nuclear regions. 



5.3. Size Distribution of H II Regions 



The size distribution of the H II regions in the original sample is displayed in Figure 
[7J showing that their sizes (diameters) range from 8 pc to 110 pc. Most of the previous 
H II region studies show that the cumulative size dist ribution of the H II regions in a 
galaxy can be well fitted with an exponential function (|van den Berghl 1 198 It iHodgd 119871 ; 



Cepa &: Beckmanlll990l ; lKnapenlll998l ; lHakobyan et al.l 120071 ). However, some other studies 



suggested that the power law form fits the differential size distribution better than the 



exponential fun c tion (IKennicutt fc Hodge 1980; E lmegreen &: Salzerlll999l ; lPleuss et al.ll2000 
Oev et alJ[2003l : lBuckalew fc Kobulnickv!l2006h . 



Figure [T^a) and (b) display the cumulative and differential size distributions of the H II 
regions, respectively. The cumulative size distribution of the H II regions with 15 pc < D < 
80 pc is fitted well with an exponential law: N(D) = N exp (—D/D Q ) where N = 77,880 
and D Q = 9.3 pc. The differential size distribution of these H II regions is fitted by a double 
power law with a break at D = 30 pc. The power law index for small H II regions with 15 
pc < D < 30 pc is «£> = —1.78 ± 0.04, whereas an = —5.04 ± 0.08 for large H II region 
with 30 pc < D < 110 pc. The small H II regions with D < 30 pc may be ionized by several 
OB stars. On the other hand, the large H II regions with D > 30 pc are considered to be 
superpositions of small H II regions, and they can be associated with stellar associations 
including several tens of OB stars. 



Figure [8]^a) shows the H II LF for small H II regions with D < 30 pc. Most of the small 
H II regions with D < 30 pc are fainter than L = 10 37 1 erg s _1 , and only 1,065 H II regions 
are brighter than L = 10 37 1 erg s -1 (solid line). Out of the H II regions with D < 30 pc, we 
selected 4,320 H II regions which are neither blended with neighbor sources nor located in 
the crowded regions (dashed line). Most of them are in the range from L = 10 35 ' 4 erg s _1 to 
10 37 2 erg s -1 . The maximum L for the isolated H II regions corresponds to about ten times 
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of M42, with an extinction-corrected value, L « 10 36 - 9 erg s _1 . On the other hand, Figure 
|8]^b) shows the H II LF for large H II regions with D > 30 pc. There are 1,808 H II regions 
brighter than L = 10 37 1 erg s _1 , and 632 H II regions fainter than L = 10 37 1 erg s _1 . 

We compared th e size distribu t ion o f the H II regions in this study with those in 
Thilker et all (bood ). IScoville et aD fj200lh and butierrez etHI fj201lh . Figure Eta) plots 
the size distributi on of the H II r egion s in the original sample (solid line), the group sample 



(dashe d line) and lScoville et al.l ( 120011 ) (dot-dashed line). In comparison with IScoville et al. 



( 120011 ). we found that there are many more small H II regions in the original sample. We 
used a single power law to fit the size distributions of H II regions for the range from 30 pc 
to 125 pc. The power law indices are a D = -5.94 ± 0.26 -3.23 ± 0.10, and -2.79 ±0.06 for 
the original sample, the group sample, and lScoville et al.l (120011 ). respectively. Thus the size 
di stribution of the orig inal sample is much steeper than that of the group sample and that 
of IScoville et al.l (120011 ) sample, while the latter two samples have similar power law indices. 



Figure Mb) shows the size distribution of t he grouped H II reg ions (solid line), the 
H II regions in the gro up sample (dashed line), in iThilker et al.l (120001 ) (dotted line) and in 
Gutierrez et al.l (120111 ) (dot-dashed line). The sizes of the H II regions in the g roup sample 
range from 8 pc to 310 pc, while those of the H II regions in lThilker et al.l ( 120001 ) range from 
80 pc to 480 pc. This large discrepancy between the two studies is considered to come from 
the different spatial resolutions. We found that t here a re many more small H II regions in the 
group sample, comparing with IGutierrez et al.l (12011). In comparison between the grouped 



H II regions and the H II regions in IGutierrez et al 



(1201 ll ). we found that although the size 



distribution for the small H II regions (D < 50 pc) have similar shapes, the size distribution 
for the large H II regions (D > 50 pc) shows different shapes, having a steeper power slope 
for the grouped H II regions. The p ower law indices are a n = —3.48 ±0.08 and —2.75 ±0.11 
for the grouped H II regions and in IGutierrez et al.l (120111 ). respectively. In other words, the 



size o f the H II regions derived from this study is generally smaller than in IGutierrez et al. 
(120 111 ). It is because we derived the sizes of the grouped H II regions fro m summing the 



numb ers of pixels contained in all individual H II regions in a group but IGutierrez et al. 



(120 111 ) considered the substantial empty volumes to measure the size of H II regions. 



Discussion 



6.1. Model Predictions for Ha Luminosity Functions 



In order to understand the nature of the observational H II LF, we performed Monte 
Carlo simulations of H II LFs, following lOey &: Clarkd (I1998I ). The Ha luminosity of an H II 
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region depends on (a) the stellar initial mass function (IMF) and (b) the relation o f Lyman 



conti nuum emission rate (QL yc ) and stellar mass. We adopted a Salpeter IMF (jSalpeter 
19551 ) with N(m*)dm* oc m*~ 2 - 35 over the range fro m 1 M Q to 100 M . We used QL yc 



depending on stellar mass given by IVacca et all ( 119961 ) which is based on stellar models with 



constraints provided by observations of individual high-mass stars. 

We considered artificial clusters with the number of stars per cluster (N*) ranging from 
1 to 10,000, and estimated the total QL yc radiating from each cluster. To reduce the statis- 
tical noise, we repeated the experiment 10,000 times, and derived the average QL yc of the 
individual cluster with a specific N*. Then we obtained the Ha luminosity of this cluster 



using L — 1.37 x 10~ 12 Q Lyc given in lScoville et all (1200 if ) 



We represented the distribution of numbers of stars in each model cluster by a power 
law: N(N*) dN* = N^dN* where N(N*)dN* is the number of stars per cluster in the range 
from AT* to iV* + dN*, and (3 is a power law index. We assumed that the star formation in 
each cluster occurs on a short time scale compared with stellar evolution time scales, and 
that all clusters in a galaxy are created at the same time. We converted simulated H II LFs 
into observational H II LFs for comparison, considering the internal extinction as adopted 
above. 

The shape of the H II LFs is determined by two parameters: and the maximum 
number of stars per cluster, N* max . To know how these two parameters affect the H II LF, 
we constructed nine H II LFs with 300,000 zero-age clusters for three power law indices 
{fi = -1.60, -2.00, and -2.40) and three values of N*, max (330, 1,000, and 10,000), plotting 
them in Figure [10J The mass of a cluster with N* jmax = 330 and (3 = 2.0 is 1,000 M and 
2,500 Mq, respectively, for the stellar mass range of 1-100 M and 0.1-100 M . The mass 
of a cluster with N* jmax = 1, 000 is 3,090 M and 7,700 M , and the mass of a cluster with 
N* >max = 10,000 is 30,900 M and 77,000 M , for the same ranges of stellar mass. The 
number of OB stars with M > 17 M Q ~ 7, ~ 20 and ~ 60 for N*, max = 330, 1, 000 and 
10,000, respectively. 

The H II LFs with N* max = 330 in the upper panels of Figure [TU] are truncated at 
L = 10 37 erg s _1 . Clusters with N* < 330 produce few H II regions brighter than L = 
10 370 erg s -1 . They show nearly flat slopes for (3 = —1.60 and -2.00, but show some 
slope for (3 = —2.40. In the middle panels, the H II LFs with N* imax = 1,000 show the 
existence of some bright H II regions with L > 10 370 erg s _1 . They show a change in slope 
at L fh 10 370 erg s _1 , which is similar to the value found for M51 H II regions in this 
study (Figure [3]). We fit the bright part (L > 10 37 ' erg s _1 ) with a power law, obtaining 
a = -1.48, -1.86, and -2.23 for (3 = -1.60, -2.00, and -2.40, respectively. Thus the values 
of a decreases as (3 decreases. The H II LFs with N* max = 10, 000 in the low panels show a 
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more obvious break at L — 10 37 ' erg s _1 . Fitting the bright part (L > 10 37 ' erg s _1 ) yields 
a = —1.55, -1.94, and -2.36 for (3 = —1.60, -2.00, and -2.40, respectively. Thus the values 
of a are very similar to those for j3, showing that the values of these two parameters get 
similar as N, increases. 

For comparison of the H II LFs in Figure [10] with the observational ones, we need to 
consider the effect of evolution of QL yc - We assumed that Ql yc remains constant during the 
main-sequence lifetimes and is zero thereaf ter. Main-sequence lifetimes (t ms ) for individual 



OB stars are obtained from iMaederl (119871 ). Figure [TT] shows the H II LFs with f3 = -2.40 
and N* jmax = 10,000 derived for four ages 0, 2, 4, and 6 Myr. We adopted a value for (3, 
similar to the observational value derived in this study. The number of bright H II regions 
decreases as ages increase. Fitting the bright part (L > 10 37,0 erg s _1 ) yields a = —2.35, 
-2.39, and -2.35 for age = 0, 2, and 4 Myr, respectively. There are few bright H II regions 
in the case of age = 6 Myr, but a fit including the faint part (L > 10 36 erg s _1 ) yields a 
similar index, a = —2.35. Thus the values of a change little depending on age, and they are 
almost the same as the input value of /3. 



6.2. A Change in Slope of Ha Luminosity Function for H II Regions and the 

Initial Cluster Mass Function 



A controversial point in the H II LF for M51 has been a break at L 1Q 38 - 8 erg s 1 shown 



in previous studies based on the ground-ba sed imaging (IKennicutt et al.l 1 19891 ; iRandl Il992 



Thilker et al.l 120001 ) . IKennicutt et al.l (119891 ) argued that the cause of this break might be a 
transition from the normal to the super giant H II regions. On the other hand , the break 
is possibly a feature caused by the low spatial resolution. IScoville et al.l (I200ll) compared 



Rand Jl992l ) based on the 



their H II LF of M51 derived from HST WFPC2 with that of 
ground-based images. They demonstrated that the high-resolution H II LF is significantly 



steeper than the low-resolution one, and showed that the break at L m 10 38 '' 
seen in their H II LF. 



erg s is not 



This break is not seen either in the H II LF for the original sample in this study, as 
shown in Figure [3j However, we can observe a break at L = io 38 - 8 erg s _1 in the H II LF 
for the group samp le in Figure |5j From the similarity between this H II LF and that in 



Thilker et al.l ( 120001 ) in Figure \MJo), we conclude that the break seen at L 
the previous studies may be due to the blending of multiple H II regions. 



10 



38.; 



erg s 



in 



On the other hand, we found a break at L — 10 erg s in the H II LF for the 
original sample as shown in Figure [3J To understand the shape of the H II LF with this 
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break we compared the observational H II LF for the original sample with the simulated 
H II LF in Figure [10J The shape of the observational H II LF is remarkably similar to 
that of the simulated H II LF in Figure [TUT i) derived for (3 = -2.40 and N* jinax = 10,000. 
The observational H II LF shows a break at L — 10 37 1 erg s -1 , which is consistent with 
the simulated H II LF with a break at L — 10 370 erg s -1 . The power law indices of the 
observational H II LF are a = —2.25 ± 0.02 for the bright part and a = —1.42 ± 0.01 
for the faint part, agreeing with those in the simulated H II LF, a = —2.36 ± 0.03 and 
a = —1.43 ± 0.02. Thus it is concluded that the break at L = 10 37,1 erg s -1 is due to 
the different N* in clusters associated with individual H II regions. In other words, the 
faint H II regions with L < 10 37,1 erg s _1 are ionized by low mass clusters with N* < 330 
(including several OB stars) or single massive stars, while the bright H II regions with 
L > 10 37 1 erg s -1 are ionized by massive clusters with N* >> 330 (including several tens 
of OB stars). 

It is noted that the cluster mass M c i is related with N*: M c \ oc N*. Therefore the 
power law index a of H II LF in the bright end is expected to be similar to the index 
(3 (-2.25) of the IMF of the clusters in M51 that may be represented by a power law: 
N(M c i) dM c i = M c i _2 ' 25 dM c i. This slope (/3) is similar to or slightly steeper than those 
derived from previous studies for M51 young clusters: from the observation of young star 
clusters in M51, the cluster m ass function is kno wn to have power la w form with its inde x 



ranging from a = -1.70±0.08 (IGieles et al.ll2006f ) to a = -2.23±0.34 flHwang fc Ledl2010h . 



6.3. Spatial Variation of Ha Luminosity Function 

Figure El shows that the shape of the H II LFs varies depending on the location within 
a galaxy. In Figure [6]^a) for the original sample, we found that although the H II LFs for the 
bright part (L > 10 37 erg s _1 ) have similar power law indices in three regions, the H II LF 
for the faint part (L < 10 37 erg s" 1 ) is steeper in the interarm region than in the arm and 
nuclear regions. In addition, there is no H II region brighter than L = 10 38 ' 2 erg s _1 in the 
interarm region. In Figure [6](b) and (c) for the grouped H II regions and the group sample, 
we also found that the H II LF is steeper in the interarm region than in the arm and nuclear 
regions. 

Although the shape of the H II LF can be affected by several factors, in cluding variation 



i n the metallicity, extinction, stellar IMF, and underlying gas dynamics, iKennicutt et al. 



(119891 ) argued that the most significant fa ctor is the ch ange in the total number of stars 



produced in typical star formation events. iRandl (119921 ) attributed the difference between 



the arm and the interarm H II LF to the existence of a more massive molecular clouds in the 
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arm region. In other words, the higher proportion of massive clouds gives rise to relatively 
more H II regions of the higher luminosity. It is reasonable to infer that the high gas surface 
densities combined with the compression by the strong densi ty wave in the arms m ay cause a 
concentration of bright H II regions on the arm region. Later lOey fc Clarkd (ll998f ) suggested 
that the interarm H II regions are on average fainter than those in the arm region because 
the ionizing sources in the arm region are younger than those in the interarm region. 

In Figure HU we show the evolution of H II LFs with f3 = -2 AO and N* max = 10,000 for 
different ages (0, 2, 4, and 6 Myr). The H II LFs shift to faintward as a function of time. 
A maximum luminosity after 4 Myr (panel c) is 0.5 dex fainter than at zero-age (panel a), 
and the power law index for the faint part (L < 10 37 ' 1 erg s _1 ) after 4 Myr is steeper than at 
zero-age. The interarm H II LF in Figure E] has a lower maximum luminosity and a steeper 
slope for the faint part than the arm H II LF. Thus the difference between the arm and the 
interarm H II LFs ca n be explained by an effect of evolution, supporting the suggestion by 



Oev k Clarkd (119981 ). 



6.4. Relation between Luminosities and Sizes of H II Regions 

Both the Ha luminosity function and the size distribution of H II regions in the original 
sample are fitted by a double power law as shown in Figure [3] and [7(b), respectively. Since 
the total Ha luminosity of an H II region is integrated from the nebular volume emission, 
it is expected that L oc D 3 . Then the power law indices (a and old) of the Ha luminosity 



funct ion and the size distribution for H II regions are related as = 2 — 3a (jOey et al. 



2003|). 



The H II LF for the original sample in Figure [3] shows that the break appears at L — 
10 371 erg s _1 , and that the power law indices are a = —1.42 ± 0.01 for the faint part 
and a = —2.25 ± 0.02 for the bright part. Using these values we derived expected values, 
old = —2.26 and old = —4.75 for the faint and bright part, respectively. On the other hand, 
the observational size distribution in Figure Ulb) shows that the power law index for the 
small H II regions with D < 30 pc is an = —1.78 ± 0.04, whereas old = —5.04 ± 0.10 for the 
large H II regions with D > 30 pc. Therefore there is a good agreement between expected 
and observational indices for the size distribution of H II regions. 

Figure H2T a) displays Ha luminosities versus sizes of the H II regions in the original 
sample. It shows a clear correlation between the two parameters. The relation between the 
luminosity and the size of H II regions is fitted by L oc £) 3 - 04±0 - 01 (solid line). This value is in 



good agreement with our expectation of L oc D 3 . It is noted, however, IScoville et all (120011 ) 
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derived a flatter relation, L oc £) 2 - 16±0 - 02 for H II regions in the central region of M51 (dashed 
line). They explained that the this flatter slope is due to blending effect: large H II regions 
are blended so that they include substantial empty or neutral volumes. T his indicates that 
this study separated blended H II regions better than IScoville et all (120011 ) . 



Figure fl2Tb) and (c) display Ha luminosities versus sizes of the grouped H II regions 
and the H II regions in the group sample, respectively. The relation between the luminosity 
and the size of H II regions for the two sa mples is fitted by L oc _D 2 - 78±0 - 02 (panel b; solid 
line) and L oc Z) 2 78±0 01 (panel c; solid line). iGutierrez et all ( 120111 ) derived a flatter relation, 
L oc D 2 - 29±om (dashed line). This discrepancy between the two studies is due to the different 
procedure used to derive the size of the H II region, as mentioned in §5.3. 



7. Summary and Conclusion 

We detected and analysed the H II regions in M51 using the wide field and high resolution 
HST ACS Ha image taken part of the Hubble Heritage program. We found about 19,60 



H II regions using the automatic H II photometry software, HIIPHOT ( IThilker et al.ll2000l ). 
Primary results are summarized as follows. 

The Ha luminosity of the H II regions ranges from L = 10 35 5 erg s _1 to 10 39 erg s" 1 . 
The H II LF is fitted by a double power law with a break at L = 10 3 erg s _1 , and the 
power law indices are a = —2.25 ± 0.02 for the bright part and a = —1.42 ± 0.01 for the 
faint part. Comparison with the simulated H II LFs suggests that this break is caused by 
the transition of H II region ionizing sources, from low mass clusters (including several OB 
stars) to massive clusters (including several tens of OB stars). 

To understand the break in the H II LF at L = 10 38 8 erg s _1 reported in the ground- 
based studies, we prod uced a group samp le, finding associations of H II regions in the original 



sample using the FOF( Da vis et al 



19851). From the comparison of the H II LF for the group 



sample with that of IThilker et al.l ( )2000l ) , we concluded that most H II regions brighter 
than L = 10 38 - 8 erg s _1 seen in the ground-based studies may be blends of multiple lower 
luminosity H II regions that are resolved as separate H II regions in this study. 

The H II LFs for the original sample have a different shape according to their positions 
in the galaxy. Although the H II LFs are fitted by a double power law with a break at 
L = 10 37 ' 1 erg s _1 , the H II LF for the faint part (L < 10 37 1 erg s _1 ) is steeper in the 
interarm region than in the arm and nuclear regions. There is no H II region brighter than 
L = 10 38 ' 2 erg s -1 in the interarm region. Observed variations in the H II LF can be explained 
by an effect of evolution of an H II region. 
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The size distribution of the H II regions in the original sample shows that the size 
of H II regions ranges from 8 pc to 110 pc in diameter. The cumulative and differential 
size distributions of the H II regions are fitted well by an exponential and power law form, 
respectively. The cumulative size distribution of the H II regions with 15 pc < D < 80 pc 
is fitted well with an exponential law: N(D) = N exp (—D/D ) with N = 77,880 and 
D = 9.3 pc. 

The differential size distribution of these H II regions is fitted by a double power law 
with a break at D = 30 pc. The power law index for the small H II regions with 15 pc < D < 
30 pc is «d = —1.78 ± 0.04, whereas «o = —5.04 ± 0.08 for the large H II region with 30 pc 
< D < 110 pc. The small H II regions with D < 30 pc may be ionized by several OB stars. 
On the other hand, the large H II regions with D > 30 pc are considered to be superpositions 
of small H II regions, and they can be associated with stellar associations involving several 
tens of OB stars. In addition, power law indices of the size distribution for the original 
sample are related with those of H II LF, and the relation between the luminosities and sizes 
of H II regions is fitted well by L oc D 3m±om . 
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Fig. 1. — A gray-scale map of the continuum-subtracted Ha image of M51. North is up and 
east to the left. The field of view is 7.2' x 10.2'. Bright discrete H II regions outline the 
spiral arms of NGC 5194, while NGC 5195 is surrounded by faint diffuse emission. 
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Fig. 2. — (a) Co mparison of H a lum inosity of 54 isolated compact H II regions common 
in this study and IScoyille et all ( 1200 ll ). The mean value of the differences is AlogL (This 
study - IScoville et al.l (120011 ))= 0.03 ± 0.12 (solid line), (b) Co mparison of Ha l umino sity of 
164 isolated compact H II regions common in this st udy and iGutierrez et al.l (120111 ). The 
mean value of the differences is AlogL (This study— IGutierrez et al.l (120111 ))= 0.01 ± 0.07 
(solid line). 
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Fig. 3. — The H II LF for the original sample (histogram). Thick dashed and thick solid 
lines represent a double power law fit for the faint part (L < 10 37 1 erg s^ 1 ) and the bright 
part (L > 10 371 erg s _1 ), respectively. The power law indices are a = —2.25 ± 0.02 for the 
bright part and a = —1.42 ± 0.01 for the faint part. An arrow marks the break point. 
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Fig. 4. — Compari son of the H II LFs fo r the original s ampl e (solid line) with those in 



Thil ker et al.l (J2000h ( dashed line in (a)), in lScoville et al.l ( 1200 ll ) (dashed line in (b)) and in 



Gutierrez et al 



(120111 ) (dashed line in (c)). Thick lines represent power law fit s . In p anel 

37.1 „.„ while Ixhilker et aD (boooh and 
much brighter than that in this 



a) and (c), this study shows a break at L = 10 erg s 



Gutierrez et al.l ( 120111 ) show a break at L — 10 38 8 erg s 



study. Arrows mark the break point. 
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Fig. 5. — (a) The H II LFs for the original (dashed line) and group (solid line) samples, (b) 
The H II LFs for the group sample (solid line), the grouped H II regions (dotted line) and in 
Thilker et al.l (120001 ) (dashed line), (c) T he H II LFs for t he gro up sample (solid line), the 
grouped H II regions (dotted line) and in iGutierrez et al.l (120111 ) (dashed line). Thick lines 
represent power law fits. Arrows mark the break point. 
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Fig. 6. — (a) The H II LFs for the arm (solid line), interarm (dotted line), and nuclear 
(dashed line) regions in the original sample, (b) The H II LFs for the arm (solid line), 
interarm (dotted line), and nuclear (dashed line) regions in the grouped H II regions, (c) 
The H II LFs for the arm (solid line), interarm (dotted line), and nuclear (dashed line) 
regions in the group sample. Thick solid lines represent power law fits. Arrows mark the 
break point. 
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Fig. 7. — The size distribution of the H II regions in the original sample, (a) Cumulative size 
distribution. The solid line represents an exponential fit. (b) Differential size distribution. 
Thick dashed and thick solid lines represent a double power law fit with a break at log 
D = 1.45 (D ~ 30 pc). The power law index for small H II regions with D < 30 pc is 
old = —5.04 ± 0.08, whereas = —1.78 ± 0.04 for large H II regions with D > 30 pc. An 
arrow marks the break point. 
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Fig. 8. — (a) The H II LFs for the small H II regions with D < 30 pc (solid line) and the 
isolated H II regions which are neither blended with neighbor sources nor located in the 
crowded regions (dashed line), (b) The H II LF for the large H II regions with D > 30 pc. 
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Fig. 9. — (a) The size distributions of the H II regi ons in the original (so lid line) and group 
(dashed line) samples in comparison with those in IScoville et al.l (120011 ) (dot-dashed line). 
Thick solid lines represent power law fits for 1.5 < log D < 2.1. (b) The size distributions 
for the grouped H II regions (soli d line) and th e H II regions in the group sample (dashed 



l ine) i n comparison with those in iThilker et al.l (120001 ) (dotted line) and in iGutierrez et al. 



(120111 ) (dot-dashed line). Thick solid lines represent power law fits for 1.7 < log D < 2.4. 



-32 - 



10000 



g 1000 



3 

Z 



100 



10 



10000 



£ 1000 



100 



10 



10000 



g 1000 



z ioo 

10 



(3 : -1.60 



IN. 



: 330 
(a) 



(3 : -1.60 




a ~ -1.48 



N, -1000 



(d) 




p : -2.00 IN, 



: 330 
(b) 



(3: -2.00 lN., raax :1000 
I (e) 



a ~ -1.86 



(3: -2.00 lN W 10000 
I (h) 




(3 : -2.40 



N», 



: 330 
(c) 



(3 : -2.40 



1000 




N» -10000 



35 36 37 38 35 36 37 38 35 36 37 38 
log L [erg s" 1 ] log L [erg s" 1 ] log L [erg s" 1 ] 



Fig. 10. — Simulated Ha LFs of the H II regions with zero age for three parent power law 
indices (j3 = -1.60, -2.00, and -2.40) and for the three maximum number of stars per cluster 
(N* jmax = 330, 1,000, and 10,000). Thick solid lines in (d) to (i) represent power law fits for 
the bright part (L > 10 37 erg s _1 ). The values of the index (a) are given in each panel. 
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Fig. 11. — Evolution of simulated H II LFs for four ages (0, 2, 4, and 6 Myr), the parent 
power law index (3 = -2.40, and the maximum number of stars per cluster N* )inax = 10,000. 
Thick solid lines represent power law fits for L > 10 37 erg s _1 . The values of the index (a) 
are given in each panel. 
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Fig. 12. — (a) The Ha luminosities ve r sus si zes of the H II regions in the original sample (dots 
and contours) and in IScoville et al.l ( 1200 ll ) (dashed line). Contours represent the number 
density of H II regions in the original sample. The solid line repre sents a power l aw fit : 
L oc £) 3 04±0 01 ) an( i the dashed line represents a fit to the data in IScoville et al.l ( 120011 ) : 
L oc £)2.i6±o.o2. (b) The Ha luminos ities versus sizes of the grouped H II regions (dots) and 
the H II regions in iGutierrez et al.l ( 1201 ll ) (dashed line). The solid line rep resents a power 



l aw fit : L oc £) 2 - 78±0 - 02 ; an d the dashed line represents a fit to the data in IGutierrez et al. 
( 1201 ll ): L oc £) 2 - 29 ± 3 , ( c ) The Ha luminos ities versus sizes of the H II regions in the group 
sample (dots) and in IGutierrez et al.l ( 1201 ll ) (dashed line). The solid line re presents a power 



l aw fit : L oc £) 2 - 78±0 - 01 ; an d the dashed line represents a fit to the data in IGutierrez et al. 
(Eml): L oc J D 2 - 2 9±o.os 
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Table 1. A catalog of M51 H II regions in the original sample 



ID 


R.A.(J2000.0) a 


Decl.(J2000.0) a 


F(Ha) b 


logL(Ha) c D d D 


Pos. c 


GID f 




[Deg] 


[Deg] 




[erg s 


H tpc] 






3 


202.3956299 


47.171803 


4.402 


36.573 


0.53 25.26 


I 




4 


202.3961639 


47.169731 


7.095 


36.780 


0.44 21.32 


I 




6 


202.3963776 


47.169849 


1.259 


36.029 


0.34 16.47 


I 




7 


202.3977661 


47.172359 


8.074 


36.836 


0.57 27.22 


I 




8 


202.3983917 


47.172512 


1.597 


36.133 


0.38 18.37 


A 




9 


202.3992462 


47.177734 


1.669 


36.152 


0.38 18.37 


A 


10 


10 


202.3993530 


47.177658 


1.290 


36.040 


0.28 13.27 


A 


10 


11 


202.3993835 


47.177444 


13.960 


37.074 


0.79 38.11 


A 


10 


13 


202.4014740 


47.182098 


13.420 


37.057 


0.58 27.75 


A 




20 


202.4026184 


47.160851 


4.296 


36.562 


0.68 32.83 


A 





a Measured in the F658N (Ha) image. 

b In the unit of [10~ 16 erg cm~ 2 s -1 ]. 

Calculated using L = 47rd 2 xF(Ha) for d=9.9 Mpc. 

d Calculated using D = 2 x (area/71-) 1 / 2 where the area is derived from the 
number of pixels contained in an H II region. 

e Position in the galaxy: A-Arm region, I-Interarm region, N-Nuclear region, 
and C-NGC 5195. 

f ID number for the group sample. 



-36 - 



Table 2. A catalog of M51 H II regions in the group sample 



ID 


R.A.(J2000.0) a 


Decl.(J2000.0) a 


F(Ha) b 


logL(Ha) c D d 


D 


Pos.° 




[Deg] 




[Deg] 




[erg s x ] 


H 


M 




3 


202.3956299 


47. 


171803 


4.402 


36.573 


0.53 


25.26 


I 


4 


202.3961639 


47. 


169731 


7.095 


36.780 


0.44 


21.32 


I 


6 


202.3963776 


47. 


169849 


1.259 


36.029 


0.34 


16.47 


I 


7 


202.3977661 


47. 


172359 


8.074 


36.836 


0.57 


27.22 


I 


8 


202.3983917 


47. 


172512 


1.597 


36.133 


0.38 


18.37 


A 


10 


202.3993530 


47. 


177654 


17.928 


37.183 


0.96 


45.96 


A 


13 


202.4014740 


47. 


182098 


13.420 


37.057 


0.58 


27.75 


A 


16 


202.4025269 


47. 


178154 


1.830 


36.192 


0.34 


16.25 


A 


17 


202.4025269 


47. 


178005 


1.471 


36.097 


0.34 


16.25 


A 


20 


202.4026184 


47. 


160851 


4.296 


36.562 


0.68 


32.83 


A 



a Measured in the F658N(Ha) image. Mean R.A.(J2000) and 
Decl.(J2000) of all individual H II regions in a group. 

b In the unit of [10 -16 erg cm~ 2 s -1 ]. Derived from summing the 
Ha fluxes of all individual H II regions in a group. 

Calculated using L = 47rd 2 xF(Ha) for d=9.9 Mpc. 

d Calculated using D = 2 x (area/71-) 1 / 2 where the area is derived from 
summing the numbers of pixels contained in all individual H II regions in 
a group. 

c Position in the galaxy: A-Arm region, I-Interarm region, N-Nuclear 
region, and C-NGC 5195. 



